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ABSTRACT 

, Context. Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) exhibit characteristic X-ray luminosi- 

ties (both soft and hard) of around 10'^'' erg s~^ and characteristic power-law, hard X-ray spectra extending to about 
200 keV. Two AXPs also exhibit pulsed radio emission. 

Aims. Assuming that AXPs and SGRs accrete matter from a fallback disk, we attempt to explain both the soft and the 
, hard X-ray emission as the result of the accretion process. We also attempt to explain their radio emission or the lack 

HH ■ of it. 

HH ' Methods. We test the hypothesis that the power-law, hard X-ray spectra are produced in the accretion flow mainly by 

(~| I bulk-motion Comptonization of soft photons emitted at the neutron star surface. Fallback disk models invoke surface 

^ ^, dipole magnetic fields of lO^'^ — 10^'^ G, which is what we assume here. 

I ■ Results. Unlike normal X-ray pulsars, for which the accretion rate is highly super-Eddington, the accretion rate is 

O ' approximately Eddington in AXPs and SGRs and thus the bulk-motion Comptonization operates efficiently. As an 

illustrative example we reproduce both the hard and the soft X-ray spectra of AXP 4U 0142-1-61 well using the XSPEC 
C/3 , package compTB. 

■ Conclusions. Our model seems to explain both the hard and the soft X-ray spectra of AXPs and SGRs, as well as their 

radio emission or the lack of it, in a natural way. It might also explain the short bursts observed in these sources. On 
1 the other hand, it cannot explain the giant X-ray outbursts observed in SGRs, which may result from the conversion 

^ , of magnetic energy in local multipole fields. 

' Key words, pulsarsdndividual (IE 1841-045, IRXS J1708-4009, 4U 0142-1-61) - X-rays: stars - stars: magnetic fields - 

0^ ■ accretion disks 

cn 
cn 

Tj" 1. Introduction ray enhancements lasting from months to more than years 

. , ,,T-, „ and are correlated, in the long term, with infrared (IR) 

O Anomalous X-ray pulsars (AXPs) and soft gamma-ray re- luminosities (Tarn et al. 2004). The transient AXPs have 

^ . peaters (SGRs) constitute a special population of young ^^^^ discovered in such X-ray outbursts (enhancements) 

' • neutron stars distmgmshed by much higher X-ray lummosi- ^^^^ ^^^^^ ^^^^^^ l^^^lg ^^^^ ^^^^^ ^^^^^^ ^^g_ 

^ ties than their rotational powers and spm periods clustered ^j^ude higher than when in quiescence (Torii et al. 1998; 

^ m a narrow range (2 - 12 s). The AXPs and SGRs are now Kouveliotou et al. 2003; Gotthelf et al. 2004; Ibrahim et 
H ; believed to belong to the same class of objects, since short 2004; Mereghetti et al. 2006, Israel et al. 2007; Muno 
. , bursts that were once believed to be a distinctive property 2007). There are both similarities and systematic dif- 
of the SGRs were observed from some of the AXPs as well ferences in the X-ray outburst light-curve morphologies of 
(Gavrul et al. 2002; Kaspi et al. 2003) . They are all spinning t.^nsient and persistent sources. Other constraints on mod- 
down with spm-period derivatives m the 10 - 10 s s ^^^^ f^^^ properties of AXPs and SGRs in the op- 
range (Woods & Thompson 2006; Kaspi 2007; Mereghetti y^^j ^^^^^ ^^^^^ relations with the X-ray 
2008 for recent reviews on AXPs and SGRs). luminosities 

In addition to these main properties, broad-band ob- 
servations have revealed many other peculiarities of these These sources are widely believed to be magnetars 
sources, which provide constraints for the models. In partic- mainly 1) because they account for the energetics of the 
ular, some of these sources are persistent with soft X-ray lu- super-Eddington soft gamma-ray bursts and 2) because 
minosities of Lx,soft ~ 10^^ - 10^^ erg s"\ while others, dis- the P and P measurements lead to large magnetic dipole 
covered in recent years, are transients with Lx,soft ~ 10^^ strengths based on the assumption that these sources 
erg s~^ in quiescence (see Table 1 of Mereghetti 2008). spin down by magnetic dipole radiation. Magnetar models 
Both transient and persistent sources show occasional X- (Duncan & Thompson 1992; Thompson & Duncan 1993; 
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1995) do explain the supcr-Eddington bursts of SGRs. 
However, the quantitative explanation of the persistent soft 
X-ray luminosity by magnetic field decay, the optical and 
IR properties during persistent states, the X-ray and the 
accompanying IR enhancements, and the period clustering 
of AXP and SGRs all seem to meet difficulties within the 
original frame of the magnetar model. 

On the other hand, fallback disk models (Chatterjee et 
al. 2000; Alpar 2001) can account for these observational 
facts. They have been developed in a series of connected, 
self-consistent ideas by means of detailed quantitative mod- 
els (Ek§i & Alpar 2003; Ertan & Alpar 2003; Ertan & Cheng 
2004; Ertan et al. 2006; Ertan & gali§kan 2006; Ertan et al. 
2007; Ertan & Erkut 2008, Ertan et al. 2009). These models 
do not explain the super-Eddington bursts. Nevertheless, 
they are compatible with the presence of magnetar fields 
provided that these fields are in higher multipoles rather 
than in the dipole component. In these models, the strength 
of the dipole magnetic field is found to be less than ~ 10^^ 
G to account for the observations. The first direct obser- 
vational support for the presence of fallback disks around 
these systems came from SPITZER observations of AXP 
4U 0142-1-61 in mid-IR bands by Wang et al. (2006). The 
same source has also been detected in optical and near-IR 
bands. Through model fits, Wang et al. (2006) show that 
the mid-IR data can be reproduced by an irradiated disk 
model. These authors propose that the ncar-IR and the 
optical luminosities have a magnetic origin, while the mid- 
IR flux originates in a passive and irradiated disk. Later 
it was shown that all the data sets from optical to mid-IR 
can be accounted for by a single disk model that is active 
and irradiated by X-rays, which themselves are produced 
by accretion of disk matter onto the neutron star (Ertan et 
al. 2007). Recently AXP IE 2259-1-586 was also detected 
in the SPITZER bands. The combined overall spectrum, 
including the earlier IR detections of this source, is similar 
to that of AXP 4U 0142+61 (Kaplan et al. 2009). 

In recent years, some of the AXPs, namely AXP IE 
1841-045 (Molkov et al. 2004; Kuiper et al. 2004), AXP 
IRXS J1708-40 (Revnivtsev et al. 2004; den Hartog et 
al. 2008a) and AXP 4U 0142-^61 (Kuiper et al. 2006; den 
Hartog et al. 2008b), have been found to emit pulsed hard 
X-rays up to 150 keV or more with isotropic luminosities 
close to the soft X-ray luminosities (see den Hartog et al. 
2008b and Mereghctti 2008). 

Detailed discussions of the magnetar models trying to 
explain the hard X-ray emission properties of AXPs can be 
found in (Heyl & Hernquist 2005; Heyl 2007; Beloborodov 
& Thompson 2007; Baring & Harding 2007). The advan- 
tages and disadvantages of the different magnetar models 
are summarized in den Hartog et al. (2008b). 

The main aim of this paper is to discuss how the hard 
and the soft X-ray components can be produced in the 
framework of accretion from a fallback disk. In doing so, 
we treat AXPs and SGRs as a class with the following rep- 
resentative properties: 

1) Luminosities are ~ 10^^ erg s~^. 

2) Soft (< 10 keV) and hard (> 10 keV) luminosities 
are similar. 

3) Dipolar magnetic field strengths are 10^^ — 10^^ G. 

4) Rotational periods are ~ 5 s. 

We find that we can explain the observations very 
naturally by considering the biilk-motion Comptonization 
(BMC) that takes place in the accretion flow above the 



polar cap. The seed photons for this Comptonization are 
provided by the polar cap, while the observed soft X-ray 
emission comes from both the polar cap and an extended 
region around it. 

As an illustrative example, we fit the observed X-ray 
spectrum (both soft and hard) of AXP 4U 0142-}-61 using 
the XSPEC package compTB (FarineUi et al. 2008). The 
fit is extremely good and the resulting parameters are very 
reasonable. 

The BMC model is presented in § 2. In § 3 we present 
the X-ray data of AXP 4U 0142-^61 and the fit to them 
using the XSPEC package compTB. In § 4 we discuss our 
model, and in § 5 we present the relation of AXPs and 
SGRs to other source classes. Finally in § 6 we give our 
conclusions. 

2. Bulk-motion Comptonization model 

2.1. Soft X-ray emission 

The soft X-ray spectra of AXPs below 10 keV can be fit- 
ted by two-component models composed of two blackbody 
spectra or a blackbody plus a steep power-law spectrum 
(Mereghetti 2008). As pointed out by Gotthelf & Halpern 
(2005), the double blackbodies are more physically moti- 
vated. The hotter component may represent photospheric 
radiation from a small hot polar cap (area ^hot), while 
the cooler component is photospheric emission from a large 
fraction of the neutron star surface (area ^cooi)- Actually, 
such a combination of spectra is also found for young ener- 
getic radio pulsars and for isolated neutron stars showing 
purely thermal emission. 

For an accreting neutron star, whose magnetic axis is 
inclined with respect to the fallback disk, the emitting re- 
gion Ahot at the base of the accretion flow will not have 
radial symmetry around the magnetic pole, but will as- 
sume the bow-shaped configuration discussed by Basko & 
Sunyaev (1976) and confirmed by the MHD calculations of 
Romanova et al. (2004; see also Bachetti et al. 2009). This 
region is heated by the infall of accreting matter and the 
produced soft X-ray photons either escape unscattered, and 
are observed as such, or get upscattered in the accretion 
fiow and produce the hard X-ray spectrum. 

2.2. Hard X-ray emission 

The equation describing upscattering of soft photons in 
a converging, ionized, fiuid fiow was first introduced and 
solved by Blandford & Payne (1981a;b) and Payne & 
Blandford (1981). Their Comptonization equation (bulk- 
motion Comptonization) was solved in the case of spheri- 
cal accretion onto a neutron star by Mastichiadis & Kylafis 
(1992), who assumed a reflective neutron-star surface. This 
last work was generalized by Titarchuk et al. (1996, 1997), 
who examined the general case of a partially reflecting in- 
ner boundary and also included a second order term in the 
flow velocity (see also Psaltis & Lamb 1997; Psaltis 2001). 

The idea in the bulk-motion Comptonization (BMC) 
model, as applied to the AXP X-ray spectra, is the follow- 
ing. Soft X-ray photons, emitted by the polar cap region, 
where the accretion occurs, find themselves in the accre- 
tion flow. For accretion rates comparable to the Eddington 
rate, the optical depth to electron scattering in the accre- 
tion flow is close to unity or more. Thus, a fraction of the 
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Nb_ = (0.93 ± 0.02) X 10^2 cm-2) from the spectral fits of 
the continuous clocking mode data. Given the uncertainties 
in the calibration of the latter mode, we consider that the 
HETG data provide a more accurate representation of the 
low-energy spectrum of 4U 0142+61. 

3.2. INTEGRAL ISGRI data 

4U 0142+61 has been extensively observed by INTEGRAL. 
We downloaded all public pointed observations from the 
INTEGRAL archive, which are longer than 3 ksec and the 
source falls within 10 degrees from the pointing direction. 
The data were analyzed with the OSA vS.cQ data analy- 
sis software provided by the Integral Science Data Analysis 
Center. We used version 8.0.1 of the instrument character- 
istics database (which provides the latest calibration data) 
and version 30.0 of the reference catalog (which provides 
a list of known sources used for the source detection and 
spectral extraction processes). First, we produced an image 
of the observed region in the 4 standard ISGRI bands (20 - 
40 keV, 40 - 60 keV, 80 - 100 keV, 100 - 200 keV) following 
the standard procedures described in the IBIS analysis user 
manuaQ. 4U 0142+61 has been clearly detected in all four 
bands. Next we extracted a spectrum of 4U 0142+61, again 
following the standard procedures for IBIS data analysis. 
As recommended, we added 2% systematic errors to the 
spectrum to account for calibration uncertainties. 

We fitted the spectrum with a power-law model, us- 
ing the latest arf and rmf files. We find an energy slope of 
r = 0.94 + 0.10, consistent with the slope ofF = 1.05 + 0.11 
reported by Kuiper et al. (2006). We also fitted the spec- 
trum with older arf files available in the instrument char- 
acteristics database in order to assess the effect of different 
calibration data on the measured spectral parameters. We 
did not find any statistically significant difference between 
the estimated spectral parameters. 

3.3. Joint Chandra and INTEGRAL ISGRI fits 

In the previous paragraphs, we have shown that our analy- 
sis of the archival Chandra and INTEGRAL ISGRI data of 
4U 0142+61 nicely reproduces the canonical model consist- 
ing of an absorbed blackbody and power-law model in the 
soft X-ray band (0.5-10.0 keV) and a hard power-law above 
20.0 keV. Next we investigated whether the same spectrum 
can be reproduced with a bulk-motion Comptonization 
model. For this reason we used the model of Farinelli et 
al. (2008) provided as an external XSPEC model (model 
COMPT53). This model includes a self-consistent treatment 
of the seed blackbody spectrum and the thermal and/or 
bulk-motion Comptonization of its photons. Therefore we 
did not include an ad-hoc blackbody component. The seed 
spectrum is described by a modified blackbody function 
S{E) oc E'^-^ /[exp{E/kTs) - 1], where kT^ is the charac- 
teristic temperature of the blackbody, and E^ is a power- 
law component that modifies the blackbody. For 7 = 3 this 
component simplifies to a pure blackbody. Since this model 
includes both bulk-motion and thermal Comptonization, 
an important parameter is the relative efficiency of the 



emitted soft photons get trapped in the flow and after sev- 
eral, nearly head-on collisions with the accreting electrons 
acquire signiflcant amounts of energy. As a result, a power- 
law, hard X-ray spectrum is produced. 

Torrejon et al. (2004) have applied the BMC model 
to the wind-accreting, slowly spinning, neutron star 4U 
2206+54, which has similar properties to our sources [L^ ^ 
2035 gj.g g-i^ 2.soft ^hard, fcJbb ~ 1.2 kcV, and a hard 
spectral tail extending up to 90 keV). The BMC model 
was also used to fit hard X-ray tails observed in neutron 
star sources by Paizis et al. (2006), and it was significantly 
improved by Farinelli et al. (2008). 

3. Application to AXP 4U 0142+61 

The AXP 4U 0142+61 is one of the brightest known AXPs. 
Therefore its spectrum can provide one of the most strin- 
gent tests of the applicability of the BMC model to the 
X-ray spectra of AXPs. It has been extensively observed 
with the Chandra, XMM-Newton, and INTEGRAL X-ray 
missions. 

To constrain the low-energy (0.5 - 10.0 keV) spectrum 
of 4U 0142+61, we opted for the available Chandra high- 
energy transmission grating (HETG) spectra. These data 
are least affected by pile-up, while they do cover the desired 
energy band. There are also Chandra data obtained in con- 
tinuous clocking mode (previously analyzed by Patel et al. 
2003); however, the calibration of this mode for spectral 
analysis is still uncertain, as is indicated by the differences 
between the published analysis of the HETG and the con- 
tinuous clocking mode data (Juett et al. 2006; Patel et al. 
2003). For the high-energy spectrum of 4U 0142+61, we use 
the extensive ISGRI INTEGRAL data, which cover the 20 
- 200 keV band. The spectral analysis has been performed 
with the XSPEC vl2.0 package (Arnaud 1996). All cited 
errors are at the 90% confidence level for one interesting 
parameter unless otherwise specified. 

3.1. Chandra data 

4U 0142+61 has been observed with the Chandra High- 
Energy Transmission Grating for 25 ksec in May 2001 
(OBSID 1018; PI C. Canizares). We used the reduced spec- 
tra available from the Chandra Grating Data Archive and 
Catalog (TGCat; http://tgcat.mit.edu/). Since even the 
0th order spectrum is att'ected by pile up, we used the 
positive and negative 1st order of the Medium Energy 
Transmission Grating and the High Energy Transmission 
Grating. The individual spectra were binned in order to 
have at least 50 counts per bin. The four spectra were fitted 
simultaneously using the relevant response matrices (rmfs) 
and ancillary responce matrices (arfs) for each order. 

A fit of the HETG data with a two component black- 
body and power-law model, affected by photoelectric ab- 
sorption by cold gas (phabs model in XPSEC) gives a 
blackbody temperature of 0.42lg'Q2 keV, a power-law slope 
of r = 3.66to'29' ^'^d an absorbing Hi column density of 
0.95lg'i[2 ^ 10^^ cm~^. These parameters are almost iden- 
tical to the parameters derived by Juet et al. (2002) from 
the analysis of the same data {kT — 0.418 ± 0.013 keV; 
r = 3.3 ± 0.4, Nu = (0.88 ± 0.13) x 10^^ cm^^). However, 
they are slightly different from the parameters derived by 
Patel et al. {kT = 0.470 ± 0.008 keV; F = 3.40 ± 0.06, 



^ http:/ /www.isdc.unige.ch/integral/download/osa_sw 
^ http:/ /isdcul3.unige.ch/Soft/download/osa/osa_doc/osa_doc- 
8.0/osa_um_ibis-7.0.pdf 
^ http:/ /heasarc. gsfc.nasa.gov/xanadu/xspec/models/comptb. html 
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Table 1. Fits of the ISGRI and Chandra data with the 
COMPTB model 



Parameter Value 

Blackbody temperature {kTs) 0.72 keV 

Index of Seed photons (7) 0.95 

Energy index of Compton spec, (a) 0.36 

Efficiency of bulk over thermal Compt. (5) 2.32 

Electron Temperature {kT^) 24.6 keV 

Ilumination factor (log(j4)) -2.3 

Normahzation (Civ)t 1-18 x 10"^ 

xVd.o.f. 166.74 / 248 



f Normalization of the seed-photon spectrum defined as 
Lsg/DiQ, where L39 is its luminosity in units of 10"^^ergs~^, 
and Dio is the distance to the source in units of 10 kpc 

two components defined a,s 6 =< i?buik >)/ < Eth >■ For 
6 = we have a pure thermal Comptonization spectrum. 
Other parameters of this model are the energy index of the 
Comptonization spectrum (a; for more details see Farinelli 
et al. 2008), the temperature of the Comptonizing electrons 
(fcTe), a factor describing the ratio between the observed 
Compton scattered spectrum and the observed seed black- 
body spectrum (A), and the normalization of the seed pho- 
ton spectrum (Cat). To model photoelectric absorption by 
cold gas, we included the XSPEC phabs model component. 

The results of this fit are presented in Table [1] and 
Figures [1] and [2] The first figure shows the data and the 
best-fit model (top panel). It also shows the ratio of the 
data to the best-fit model (bottom panel) . Figure [2] shows 
the model and the data corrected for instrumental effects 
(unfolded spectrum) in E^f{E) space. It is evident from 
this figure that the model gives a very good fit to the 
broad band spectrum of 4U 0142-f61 (xl = 166.7/248). 
The temperature of the seed photons {kT = 0.72 keV) is 
slightly higher than the temperature estimated from the 
blackbody fit of the Chandra data alone. However, the pa- 
rameter 7 = 0.95, which modifies the seed blackbody spec- 
trum indicates that it is not a pure blackbody. Furthermore, 
we find that the bulk-motion Comptonization dominates 
over thermal Comptonization (6 = 2.32) and that only 
0.3% of the observed photons have been Compton scattered 
[log(^) = —2.3]. These results show that the spectrum of 
4U 0142-1-61 can be reproduced equally well with a self- 
consistent BMC model with very reasonable parameters. 

4. Discussion 

The advantage of the bulk-motion Comptonization model 
is that it uses a proven concept, which has been applied 
to many cases including the slowly pulsating source 4U 
2206-f54, which in turn shows a low-luminosity (10'^^ erg 
s"^) and a hard X-ray tail like the AXPs (Torrejon et 
al. 2004). Another advantage is that the whole spectrum, 
from the soft thermal component to the hard X-ray tail, 
is explained by a single, well-understood mechanism us- 
ing as few parameters as possible. The accreting matter 
both heats the polar cap, which emits the seed photons for 
Comptonization, and provides the high-energy electrons, 
which do the Comptonization. Depending on detailed con- 
ditions, both the soft and the hard X-ray components have 
luminosities that are within the same order of magnitude 
as observed. Since the soft and hard components come from 
the same place at the neutron star surface, their maxima 
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Fig. 1. Top panel: The Chandra HEG and MEG (first pos- 
itive order shown with red and green points respectively) 
and INTEGRAL ISGRI data (shown in black) along with 
the best-fit COMPTB model. Bottom panel: Ratio between 
the data and the best-fit model. 
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Fig. 2. The best-fit COMPTB model, along with the data 
corrected for instrumental effects shown in E'^f{E) space. 



occur at basically the same rotational phase. This agrees 
well with the observational fact that these sources show a 
very slow change in pulse profile with energy (Kaspi 2007). 
On the other hand, there are fine structures in the energy- 
dependent pulse profiles, which can be explained by cy- 
clotron resonance efi'ects, as in the case of normal X-ray 
pulsars having dipole magnetic fields of similar strengths 
to the ones postulated by our model for AXPs and SGRs. 

The radio emission observed in two of the AXPs [IE 
1547-54 (Camilo et al. 2007) and XTE J1810-197 (Halpern 
et al. 2005; Camilo et al. 2006)] can also be explained by 
the BMC model. Usually, it is thought that the radio emis- 
sion of a neutron star is quenched by the accretion, but 
there may be little or no overlap in our model between the 
circumpolar zone of radio emission and the bow-shaped ac- 
cretion zone (Romanova et al. 2004). Indeed, according to 
the canonical expression, the size of the radio-active polar 
cap is Tpc = 10'*P~^/^ ^ 5 X 10'^ cm. It seems likely that 
the bow-shaped accretion region is farther from the mag- 
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netic pole. Therefore, both radio and X-ray emission may 
coexist. 

A necessary requirement for a detectable radio flux 
is that the magnetic dipole power E is strong enough. 
Comparing E = {l/6c^)B'^R'^{2Tr/P)^, where B is the sur- 
face magnetic field, R the radius of the neutron star, P the 
rotational period and c the speed of light, with the observed 
radio luminosities of the two pulsars, we estimate the ra- 
dio emission efficiencies to be ~ 10^"* and ^ 0.03 for IE 
1547-54 and XTE J1810-197, respectively, when assuming 
a polar magnetic field of 10^^ G. A lower magnetic field 
strength is also allowed in our model, because the observed 
radio luminosity is not isotropic. These efficiencies are quite 
large, but comparable to those of the most efficient normal 
radio pulsars with periods > 2 s, taken from the ATNF 
catalog http: / / www. atnf . csiro . au / research / pulsar / psrcat /, 
(Manchester et al. 2005). 

The fact that most AXPs and SGRs have remained un- 
detected in the radio band despite deep searches (Burgay 
et al. 2007) may be due to their long periods, which may 
place them beyond the death line. Another reason could 
be that long-period radio pulsars have narrow beams with 
widths scaling with P~^^^ (Lyne & Manchester 1988). For 
a period of 6 s, the full beam is then only ~ 7 degrees. 
Therefore the probability of detecting the radio beam of an 
AXP is only a few percent. 

5. Relation to other source classes 

The main point of this paper is to explain the hard X- 
ray spectra of AXPs and SGRs as a result of accretion 
from a fallback disk. In addition, we attribute the soft X- 
ray spectra of AXPs and SGRs to thermal, photospheric 
radiation from the polar cap (hotter component) and from 
the bulk surface of the neutron star (cooler component). 
This latter point leads us to the following thoughts: 

The class of isolated neutron stars showing purely ther- 
mal emission that was discovered by ROSAT (XDINS, 
Magnificent Seven) exhibits blackbody emission as well, 
though at lower temperatures (by a factor ~ 10) and lower 
luminosities (by ^ 10^) compared with AXPs. Like the 
latter, the XDINs have long periods and strong magnetic 
fields (> 10^3 G, Haberl 2007). Their period derivatives 
and proper motion data indicate ages of a few hundred 
thousand years, and they are not associated with super- 
nova remnants. All these facts show that XDINs are older 
than AXPs, supporting the idea that some (or possibly all) 
XDINs have been AXPs before they exhausted their resid- 
ual disk and have subsequently cooled down. 

Another young neutron-star population, namely the 
central compact objects (CCOs, see review by Pavlov et 
al. 2004), are probably related to AXPs, SGRs, and XDINs 
as well. It was originally proposed by Alpar (2001) that sim- 
ilarities and differences in these neutron star classes could 
be explained if fallback disks with different properties are 
included in the initial parameters of the evolution models 
in addition to initial spin period and the magnetic moment 
of the neutron star. 

6. Conclusions 

Contrary to common belief, AXPs and SGRs may not be 
that different from normal accreting X-ray pulsars. The ma- 



jor differences between the two classes of X-ray sources are 
the extent of the hard X-ray spectrum and the X-ray lu- 
minosity, which is lower for AXPs/SGRs by two to three 
orders of magnitude. The hard X-ray power-law spectral 
index is essentially the same for both classes, but the spec- 
trum of normal X-ray pulsars extends to about 20 keV, 
while that of the anomalous X-ray pulsars extends to about 
200 keV. In our model, this difference is caused by the nor- 
malized (to the Eddington value) accretion rate being near 
unity in AXPs/SGRs, and bulk- motion Comptonization 
takes place and produces the hard X-ray tail, which ex- 
tends to ~ 200 keV. 

Our model naturally explains the ^ 100 % pulsations 
observed at ~ 100 keV, because they come from the accre- 
tion flow above the polar cap. It also explains the similarity 
between AXPs/SGRs and XDINS, which has often been 
noted. It may be that XDINS are evolved AXPs with no 
remaining accretion disk. The absence of cyclotron lines in 
the observed X-ray spectra is a natural consequence of our 
model. For a dipole magnetic field of ~ 10^^ G, the elec- 
tron cyclotron line would appear a.t E > 100 keV, where 
the photon statistics are not good enough. On the other 
hand, a proton cyclotron line would appear at i? ~ 0.05 
keV, where interstellar absorption is huge. 

Our model could in principle explain the short bursts 
observed from AXPs and SGRs as the result of infalling 
lumps of matter with appropriate mass. We defer an anal- 
ysis of this process to another paper. On the other hand, 
our model does not explain the giant bursts observed from 
AXPs and SGRs. These may indeed require very strong, 
magnetar-type magnetic fields. Our model does not exclude 
the existence of such super strong magnetic fields if they re- 
side in multipole components. 

Last but not least, we feel that the detection of a fall- 
back disk around the neutron stars in AXP 4U 0142-1-61 
and AXP IE 2259-1-586 gives strong observational support 
to our model. 
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